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Abstract 

We show that axion domain walls gain spontaneous magnetization in early 
universe by trapping either electrons or positrons with their spins polarized. 
The reason is that the walls produces an attractive potential for these parti- 
cles. We argue that the wall bounded by an axionic superconducting string 
leaves a magnetic field after its decay. We obtain a field ~ 10~ 23 Gauss on 

the scale of horizon at the recombination. 
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Magnetic fields in galaxies or stars are present in our universe. Nevertheless, their origin is 
still unknown although there are several cosmological origins proposed ; it could arise 
during quark-hadron phase transition electro- weak phase transition |2],|3]], inflation Q 
e.t.c. It is generally assumed that dynamo process |]] amplifies the seed of small primordial 
magnetic field ~ 10~ 18 Gauss generated by these origins to the observed one 10~ 6 Gauss in 
galaxy. 

Previously we have discussed |§ a magnetic field associated with domain walls, which gain 
ferromagnetism owing to fermion zero modes bounded to the walls. But the zero modes 
does not necessarily exist in any realistic models. In this letter we show that axion domain 
walls, although they do not possess such zero modes, become ferromagnetic by trapping 
either electrons or positrons in the hot universe. The walls leave dipole magnetic fields after 
their decay with various strengths and sizes, e.g. a field of ~ 10 20 Gauss on scales of ~ 10~ 3 
cm at the temperature ~ 100 MeV. These dipole fields generate a large scale magnetic field. 
The field on the scale of horizon is estimated to possess its strength 10 -23 Gauss at the 
recombination of photons and electrons. This magnetic field is a candidate of primordial 
magnetic fields which lead to galactic or intergalactic magnetic field at present. 

Let us begin to briefly explain axion domain walls H. The axion is a Goldstone mode H] 



of Pecci-Quinn global U(l) symmetry |]T(| which is broken spontaneously with the energy 
scale fpQ, 10 10 ~ 10 12 GeV. This massless axion gains a mass m a through effects of QCD 
instantons due to anomaly [|ll[] of the symmetry. Since the axion field a is essentially the 
phase 9 of a Higgs field a ~ fpQ& %6 ( a = fpq9 ), the vacuum 9 = of the axion is degen- 
erate with the vacuum 9 = 2tt. This fact leads to domain walls between these degenerate 
vacua. Namely the domain walls are produced during QCD phase transition at which QCD 
instantons work effectively. Besides the domain wall solitons, string solitons associated with 
the breaking of the U(l) symmetry are produced at the temperature Jpq. In general they 
are superconducting ||12| . Since we consider an axion model with the color anomaly of the 
Pecci-Quinn symmetry equal to one, domain walls are surrounded by these strings. Fur- 
thermore, they decay soon after their production without dominating the universe. As we 



show below, however, the walls leave magnetic fields with various sizes and strengths. This 
is because the walls gain magnetic moments by trapping either electrons or positrons with 
their spins polarized. It turns out that in the decay of the walls the superconducting strings 
play important roles on producing the magnetic fields. ( For definiteness we use numerical 
values such that f PQ = 10 12 GeV and m a = 10~ 5 eV throughout the paper. ) 

In order to see the spontaneous magnetization of the wall we first show that the axion 
domain wall generates an attractive potential for electrons and positrons with a particular 
polarization. Suppose an axial vector coupling between electron field tp and axion field a, 

Lm = gd a a^^i) (1) 

where g is the coupling constant whose value is given by cfpq, with a numerical positive 

constant c; the value of c depends on models of the axion. Then, we can see easily that the 

wall produce a spin dependent potential for electrons and positrons. Assuming that the flat 

wall is located at x 3 = 0, we rewrite Dirac equation of ip — (^j such that 

d 2 ( d 2 \ 

£U = ^~2m + gd ^ aa ^> u [ £v = (~2m~ ~ gd^aa^v 1 (2) 

with energy E = m + e(E = — m — e ) ; m is the mass of electron. We have taken the 

nonrelativistic limit so that only large component u ( v ) is exhibited; small component v ( 

u ) is given by v — —id 3 u/2m ( u = id 3 v/2m ). 

We note that the typical scale of spatial variation of 6 is given by the axion mass m a ; it 
increases from to 2n along x 3 axis. Thus it turns out that the potential gd 3 aa 3 for both 
of electrons and positrons is attractive for states with spin down, and repulsive for states 
with spin up. It has approximately the width of m^ 1 ~ 1 cm and the depth of m a ~ 1CT 5 
eV when the constant of c is order of one. Hence the potential depth is quite shallow but 
the width is so large that the potential can accommodate bound states for electrons and 
positrons. Obviously, spins of these particles bounded to the wall are aligned. This causes 
the ferromagnetism of the axion domain wall. 

It seems that in the universe with its temperature ~ 100 MeV the states with such 
small binding energies are irrelevant for the property of the wall. But we should note that 



there exists a small fraction of the particles occupying the bound states even in such hot 
universe; there is tiny but nonvanishing probability of a particle occupying the state. This 
small fraction of the particles gives rise to the magnetic property of the wall whose relevant 
energy scale is much smaller than the temperature. 

Secondly we show that for a strong magnetic field either electrons or positrons have no 
bound states even in such attractive potential. For the purpose let us see the energy spectra 
of these particles with spins down under the magnetic field imposed perpendicular to the 
wall. When the field is pointed to the direction along positive x 3 axis, the binding energies 
of these particles are found easily, 

£ ^_ m „ + ii +w „ forelectron 

2m 
k 2 

£k,n+i = _ m a + 7p- + w(n + 1) for positron (3) 

Zt lit 

where n is integer ( > ) and uj is the cyclotron frequency ( uj = eB/m ); k 3 is momentum 
along rr 3 axis. For convenience we have simplified the attractive potential by a square well 
potential with a width equal to m^ 1 and a depth equal to m a . The difference in the energy 
spectra between electrons and positrons comes from the difference in the directions of the 
magnetic moments of these particles in the attractive potential of the wall. 

We note that when binding energies are positive ( e > ), the particles are not bounded. 
Such localized states are quite unstable against any small perturbations, e.g. couplings with 
scattering states of the particles, oscillation modes of the wall e.t.c. Especially thermal 
fluctuations would destroy such unstable states. Thus the states with e > must decay. 
Hence for the sufficiently strong magnetic field eB > m a m, positrons do not form bound 
states on the wall since the lowest energy Sk=o,\ ( = —fn a + u; ) of positrons becomes 
positive. Similarly, when the magnetic field is pointed to the direction along the negative 
x 3 axis, electrons do not form bound states. 

If the magnetic field is generated by these polarized electrons or positrons themselves 
bounded to the wall, it means that the wall gains magnetization spontaneously. As will 
be shown, when the temperature is higher than 180 MeV, the wall gains the spontaneous 
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magnetization with sufficiently strong magnetic field. 

Now we show the spontaneous magnetization of the wall by calculating free energy of the 
gas of electrons bounded to the wall. The free energy Fm of the system under the magnetic 
field is defined as 

F M = -(3~ 1 ]T log (1 + e -( m +^) (4) 

£fe,n<0 

where we have taken account of only contribution of electrons bounded to the wall, since 
positrons have no bound states under the sufficiently strong magnetic field. Summation, 
J2 £kn <o implies N d J2 n =o I m a l dkz/ (2ir) with the condition that e^ n < 0; N d is the degener- 
acy of each Landau level given by eBS/2ir where S is the surface area of the wall, fi^ 1 is 
the temperature of the universe. We have assumed the chemical potential to be vanishingly 
small. Noting that only states relevant in the summation are states of the lowest Landau 
level, we find that F M is approximately given by 



m- l N d 2\og2 ^2^T a eBS\og2 yJ2m/m a 
Fm ~ 2^ " 2^ (5) 

in the limit of the temperature being much higher than the mass of electron, » m. 
Thus magnetization of the wall is 



m a dF M _ e log 2 yJ2mm a 

S dB 2n 2 [3 1 ' 

This magnetization induces a current J at the boundary of the wall; J = M/m a . Then, the 
current induces a magnetic field which is not uniform; but it is approximated to B = J/R 
where R is the radius of the wall. We identify this magnetic field B with the field B in the 
free energy Fm- In this way we obtain the magnetic field generated by polarized electrons 
bounded to the wall, 



_ ey/21og2 Vm 

For the consistency this magnetic field has to satisfy a condition, eB > mm a . Otherwise 
the assumption that only particles bounded to the wall are electrons and that positrons have 



no contribution to the free energy, does not hold. Therefore when the condition is satisfied, 
the wall may gain the spontaneous magnetization by trapping either electrons or positrons. 
As have been shown in numerical simulations |H| , energetically important domain walls 



are ones with boundaries of horizon size. Thus we equate R with the distance to the 
horizon. Noting that R is about equal to 0.3M pi /? 2 / y/J ( M v \ is Planck mass and / is 
dynamical massless degrees of freedom at temperature ) |T^| , we find from eq(^) that 
the consistency condition, 

B > mm a /e ~ 230Gauss (8) 

is satisfied when the temperature is higher than 180 MeV. Therefore the wall gain the 
magnetization M spontaneously when it is produced during QCD phase transition ( ~ 200 
MeV ). Note that m a in the above formulae represents the depth of the potential when the 
coupling constant g in eq(|l]) between the axion and electron is the order of fpg. On the 
other hand if g is the order of ag ( a = e 2 /47r ) as predicted in a model of a hadronic axion 
PI, then the depth of the potential m a is replaced by am a . Thus the minimum temperature 
needed for the appearance of the magnetization is given approximately by 15 MeV. 

In order to establish the fact that the spontaneous magnetization actually arises, we 
have to make sure that the above free energy Fm is lower than that of the state with no 
magnetization. Such a free energy is given by 

F = -2/T* E log(l + e~^) ~ -2^1og2 £ = _ ^o g 2™rnf^ a 

with Ek = —m a + k 2 /2m, where we have taken both contributions from electrons and 
positrons, and have taken the limit of the temperature being much larger than m; k is 
3-momentum. Comparing this free energy F with the above one Fm in eq(|5|) we find that 
Fq > Fm when a condition, eB > Amm a /3, is satisfied; the condition is roughly the same as 
the above condition eq(|^). Therefore spontaneous magnetization of the axion domain wall 
of horizon size arises when the temperature is higher than about 180 MeV. In the case the 
magnetic field stronger than 200 Gauss is generated. 



It seems apparently strange that the spontaneous magnetization arises in the side of 
higher temperature than the critical one. In general, ordered states appear in the side of 
lower temperature than the critical one, possessing smaller entropies than those of disordered 
states. In our case, however, the entropy of the state with the magnetization is larger than 
that of the state without the magnetization, when B is larger than a critical value given 
in eq(||). This is because the number of the bound states, especially the states of the 
lowest Landau level increase with the magnetic field ( the degeneracy Na increases with 
B ) and consequently the entropy becomes large. This magnetic field is proportional to 
the magnetization which increases with the temperature: The thermal fluctuation leads to 
large orbital angular momentum and hence large magnetic moment. ( Note that states with 
various angular momenta are degenerate in a Landau level. ) Thus when the temperature 
becomes large, the magnetic field increases and the entropy becomes large. Eventually it 
dominates over the entropy of the state without the magnetization and consequently the 
spontaneous magnetization can arise. 

As we have shown, the axion domain walls trap either electrons or positrons with their 
spins polarized and gain the magnetization. Although they also gain electric charges, the 
charges are screened immediately due to large electric conductivity of the universe. Thus 
charge neutrality is kept. 

Finally we discuss magnetic fields left after the decay of the walls. We easily understand 
that during the decay of the walls, the magnetic field becomes strong as the radius of the 
walls decreases. To see it we note that the walls are surrounded by superconducting strings 
which must carries the boundary current J = M/m a . This current is supposed to be carried 
by fermion zero modes on the strings. Then as the number of the zero modes ( ~ JR 

) is conserved, the current increases as the radius R of the string decreases. The strings 



shrink with increasing the current until the current is saturated |I2| with its maximal value 
qruf/Tr; q and rrif are the electric charge and the mass of the fermion making zero modes on 
the strings. Subsequently, the strings shrink without increasing the current but by emitting 
the fermion zero modes. Therefore, when we start with, for example, a magnetic field 10 2 



Gauss of the walls on the scale of horizon Rh ~ 10 6 cm at the temperature 100 MeV, the field 
B c = 10 2 (Rh/L c ) 2 ~ 10 20 Gauss on the scale of L c is achieved when the current is saturated 
with the radius of the strings being equal to L c . L c is determined by MRh/m a = qmfL c /7T 
( conservation of the zero modes ); L c ~ 10~ 3 cm with use of m/ = 10 12 GeV and q — 1 
( this is a typical energy scale Jpq of the fermion mass which is generated through Higgs 
mechanism associated with Pecci-Quinn symmetry ). For smaller radius R of the strings, 
stronger magnetic fields ( oc 1/R ) arise on smaller scales of R. On the other hand the 
large scale magnetic field should be determined by these randomly oriented dipole fields. 
Thus they produce a magnetic field B(L) on the scale of L such that B(L) = B C (L C / L) 3//2 
Hill. The field evolves to the field B re (L re ) = S(L)(leV/100MeV) 2 on the scale of 
L re = L(100MeV/leV) at the recombination ( 1 eV ). Note that the magnetic flux is 
conserved in the early universe owing to large electric conductivity; Ba 2 RW = constant ( clrw 
is the cosmic scale factor in Robertson- Walker metric ). As we are concerned with the field 
on scales of horizon size ~ O.lMpc at the recombination, we find that B re (Lh) = 10~ 23 
Gauss. Similarly we obtain the magnetic field ~ 10 8 Gauss on scales of ~ 10 4 cm at the 
nucleosynthesis. This is sufficiently small not to affect seriously the production of light 
elements [fL7| . In this way the ferromagnetic axion domain walls generates the primordial 
magnetic fields. 
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